ChIATAC is an efficient strategy for multi-omics mapping of 3D epigenomes from low-cell inputs.

Connecting genes to their cis-regulatory elements has been enabled by genome-wide mapping of chromatin interactions using proximity ligation in ChIA-PET, Hi-C, and their derivatives. However, these methods require millions of input cells for high-quality data and thus are unsuitable for many studies when only limited cells are available. Conversely, epigenomic profiling via transposase digestion in ATAC-seq requires only hundreds to thousands of cells to robustly map open chromatin associated with transcription activity, but it cannot directly connect active genes to their distal enhancers. Here, we combine proximity ligation in ChIA-PET and transposase accessibility in ATAC-seq into ChIATAC to efficiently map interactions between open chromatin loci in low numbers of input cells. We validate ChIATAC in Drosophila cells and optimize it for mapping 3D epigenomes in human cells robustly. Applying ChIATAC to primary human T cells, we reveal mechanisms that topologically regulate transcriptional programs during T cell activation.

Oncogenic extrachromosomal DNA functions as mobile enhancers to globally amplify chromosomal transcription.

Extrachromosomal, circular DNA (ecDNA) is emerging as a prevalent yet less characterized oncogenic alteration in cancer genomes. We leverage ChIA-PET and ChIA-Drop chromatin interaction assays to characterize genome-wide ecDNA-mediated chromatin contacts that impact transcriptional programs in cancers. ecDNAs in glioblastoma patient-derived neurosphere and prostate cancer cell cultures are marked by widespread intra-ecDNA and genome-wide chromosomal interactions. ecDNA-chromatin contact foci are characterized by broad and high-level H3K27ac signals converging predominantly on chromosomal genes of increased expression levels. Prostate cancer cells harboring synthetic ecDNA circles composed of characterized enhancers result in the genome-wide activation of chromosomal gene transcription. Deciphering the chromosomal targets of ecDNAs at single-molecule resolution reveals an association with actively expressed oncogenes spatially clustered within ecDNA-directed interaction networks. Our results suggest that ecDNA can function as mobile transcriptional enhancers to promote tumor progression and manifest a potential synthetic aneuploidy mechanism of transcription control in cancer.

Reduced subgenomic RNA expression is a molecular indicator of asymptomatic SARS-CoV-2 infection.

Background: It is estimated that up to 80% of infections caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are asymptomatic and asymptomatic patients can still effectively transmit the virus and cause disease. While much of the effort has been placed on decoding single nucleotide variation in SARS-CoV-2 genomes, considerably less is known about their transcript variation and any correlation with clinical severity in human hosts, as defined here by the presence or absence of symptoms. Methods: To assess viral genomic signatures of disease severity, we conducted a systematic characterization of SARS-CoV-2 transcripts and genetic variants in 81 clinical specimens collected from symptomatic and asymptomatic individuals using multi-scale transcriptomic analyses including amplicon-seq, short-read metatranscriptome and long-read Iso-seq. Results: Here we show a highly coordinated and consistent pattern of sgRNA expression from individuals with robust SARS-CoV-2 symptomatic infection and their expression is significantly repressed in the asymptomatic infections. We also observe widespread inter- and intra-patient variants in viral RNAs, known as quasispecies frequently found in many RNA viruses. We identify unique sets of deletions preferentially found primarily in symptomatic individuals, with many likely to confer changes in SARS-CoV-2 virulence and host responses. Moreover, these frequently occurring structural variants in SARS-CoV-2 genomes serve as a mechanism to further induce SARS-CoV-2 proteome complexity. Conclusions: Our results indicate that differential sgRNA expression and structural mutational burden are highly correlated with the clinical severity of SARS-CoV-2 infection. Longitudinally monitoring sgRNA expression and structural diversity could further guide treatment responses, testing strategies, and vaccine development.

Chromatin interaction analyses elucidate the roles of PRC2-bound silencers in mouse development.

Lineage-specific gene expression is modulated by a balance between transcriptional activation and repression during animal development. Knowledge about enhancer-centered transcriptional activation has advanced considerably, but silencers and their roles in normal development remain poorly understood. Here, we performed chromatin interaction analyses of Polycomb repressive complex 2 (PRC2), a key inducer of transcriptional gene silencing, to uncover silencers, their molecular identity and associated chromatin connectivity. Systematic analysis of cis-regulatory silencer elements reveals their chromatin features and gene-targeting specificity. Deletion of certain PRC2-bound silencers in mice results in transcriptional derepression of their interacting genes and pleiotropic developmental phenotypes, including embryonic lethality. While some PRC2-bound elements function as silencers in pluripotent cells, they can transition into active tissue-specific enhancers during development, highlighting their regulatory versatility. Our study characterizes the molecular profile of silencers and their associated chromatin architectures, and suggests the possibility of targeted reactivation of epigenetically silenced genes.

Ultra-long Read Sequencing for Whole Genomic DNA Analysis.

Third generation single-molecule DNA sequencing technologies offer significantly longer read length that can facilitate the assembly of complex genomes and analysis of complex structural variants. Nanopore platforms perform single-molecule sequencing by directly measuring the current changes mediated by DNA passage through the pores and can generate hundreds of kilobase (kb) reads with minimal capital cost. This platform has been adopted by many researchers for a variety of applications. Achieving longer sequencing read lengths is the most critical factor to leverage the value of nanopore sequencing platforms. To generate ultra-long reads, special consideration is required to avoid DNA breakages and gain efficiency to generate productive sequencing templates. Here, we provide the detailed protocol of ultra-long DNA sequencing including high molecular weight (HMW) DNA extraction from fresh or frozen cells, library construction by mechanical shearing or transposase fragmentation, and sequencing on a nanopore device. From 20-25 µg of HMW DNA, the method can achieve N50 read length of 50-70 kb with mechanical shearing and N50 of 90-100 kb read length with transposase mediated fragmentation. The protocol can be applied to DNA extracted from mammalian cells to perform whole genome sequencing for the detection of structural variants and genome assembly. Additional improvements on the DNA extraction and enzymatic reactions will further increase the read length and expand its utility.

Complex Microbiome in Brain Abscess Revealed by Whole-Genome Culture-Independent and Culture-Based Sequencing.

Brain abscess is a severe infectious disease with high mortality and mobility. Although culture-based techniques have been widely used for the investigation of microbial composition of brain abscess, these approaches are inherent biased. Recent studies using 16S ribosomal sequencing approaches revealed high complexity of the bacterial community involved in brain abscess but fail to detect fungal and viral composition. In the study, both culture-independent nanopore metagenomic sequencing and culture-based whole-genome sequencing using both the Illumina and the Nanopore platforms were conducted to investigate the microbial composition and genomic characterization in brain abscess. Culture-independent metagenomic sequencing revealed not only a larger taxonomic diversity of bacteria but also the presence of fungi and virus communities. The culture-based whole-genome sequencing identified a novel species in Prevotella and reconstructs a Streptococcus constellatus with a high GC-skew genome. Antibiotic-resistance genes CfxA and ErmF associated with resistance to penicillin and clindamycin were also identified in culture-based and culture-free sequencing. This study implies current understanding of brain abscess need to consider the broader diversity of microorganisms.

Mapping the Global Chromatin Connectivity Network for Sox2 Function in Neural Stem Cell Maintenance.

The SOX2 transcription factor is critical for neural stem cell (NSC) maintenance and brain development. Through chromatin immunoprecipitation (ChIP) and chromatin interaction analysis (ChIA-PET), we determined genome-wide SOX2-bound regions and Pol II-mediated long-range chromatin interactions in brain-derived NSCs. SOX2-bound DNA was highly enriched in distal chromatin regions interacting with promoters and carrying epigenetic enhancer marks. Sox2 deletion caused widespread reduction of Pol II-mediated long-range interactions and decreased gene expression. Genes showing reduced expression in Sox2-deleted cells were significantly enriched in interactions between promoters and SOX2-bound distal enhancers. Expression of one such gene, Suppressor of Cytokine Signaling 3 (Socs3), rescued the self-renewal defect of Sox2-ablated NSCs. Our work identifies SOX2 as a major regulator of gene expression through connections to the enhancer network in NSCs. Through the definition of such a connectivity network, our study shows the way to the identification of genes and enhancers involved in NSC maintenance and neurodevelopmental disorders.

Multiplex chromatin interactions with single-molecule precision.

The genomes of multicellular organisms are extensively folded into 3D chromosome territories within the nucleus1. Advanced 3D genome-mapping methods that combine proximity ligation and high-throughput sequencing (such as chromosome conformation capture, Hi-C)2, and chromatin immunoprecipitation techniques (such as chromatin interaction analysis by paired-end tag sequencing, ChIA-PET)3, have revealed topologically associating domains4 with frequent chromatin contacts, and have identified chromatin loops mediated by specific protein factors for insulation and regulation of transcription5-7. However, these methods rely on pairwise proximity ligation and reflect population-level views, and thus cannot reveal the detailed nature of chromatin interactions. Although single-cell Hi-C8 potentially overcomes this issue, this method may be limited by the sparsity of data that is inherent to current single-cell assays. Recent advances in microfluidics have opened opportunities for droplet-based genomic analysis9 but this approach has not yet been adapted for chromatin interaction analysis. Here we describe a strategy for multiplex chromatin-interaction analysis via droplet-based and barcode-linked sequencing, which we name ChIA-Drop. We demonstrate the robustness of ChIA-Drop in capturing complex chromatin interactions with single-molecule precision, which has not been possible using methods based on population-level pairwise contacts. By applying ChIA-Drop to Drosophila cells, we show that chromatin topological structures predominantly consist of multiplex chromatin interactions with high heterogeneity; ChIA-Drop also reveals promoter-centred multivalent interactions, which provide topological insights into transcription.

Linked-read Sequencing Analysis Reveals Tumor-specific Genome Variation Landscapes in Neurofibromatosis Type 2 (NF2) Patients.

HYPOTHESIS: We hypothesize that genomic variants including deletions, insertions, inversions, and tandem duplications beyond the changes in tumor suppressor NF2 gene affect gene expression of tumor-specific pathways in vestibular schwannomas (VS) patients with Neurofibromatosis type 2 (NF2), thus contributing to their clinical behavior. BACKGROUND: Genomic variation could reconfigure transcription in NF2 transformation process. Therefore, genome-wide high-resolution characterization of structural variants (SV) landscapes in NF2 tumors can expand our understanding of the genes regulating the clinical phenotypes in NF2-associated VS. METHODS: We performed whole-genome haplotype-specific structural variation analysis using synthetic linked reads generated through microfluidics-based barcoding of high molecular weight DNA followed by high-coverage Illumina paired-end whole-genome sequencing from 10 patients' tumors of different growth rates and their matching blood samples. RESULTS: NF2 tumor-specific deletions and large SVs were detected and can be classified based on their association with tumor growth rates. Through detailed annotation of these mutations, we uncover common alleles affected by these deletions and large SVs that can be associated with signaling pathways implicated in cell proliferation and tumorigenesis. CONCLUSION: The genomic variation landscape of NF2-related VS was investigated through whole-genome linked-read sequencing. Large SVs, in addition to deletions, were identified and may serve as modulators of clinical behavior.

Picky comprehensively detects high-resolution structural variants in nanopore long reads.

Acquired genomic structural variants (SVs) are major hallmarks of cancer genomes, but they are challenging to reconstruct from short-read sequencing data. Here we exploited the long reads of the nanopore platform using our customized pipeline, Picky ( https://github.com/TheJacksonLaboratory/Picky ), to reveal SVs of diverse architecture in a breast cancer model. We identified the full spectrum of SVs with superior specificity and sensitivity relative to short-read analyses, and uncovered repetitive DNA as the major source of variation. Examination of genome-wide breakpoints at nucleotide resolution uncovered micro-insertions as the common structural features associated with SVs. Breakpoint density across the genome is associated with the propensity for interchromosomal connectivity and was found to be enriched in promoters and transcribed regions of the genome. Furthermore, we observed an over-representation of reciprocal translocations from chromosomal double-crossovers through phased SVs. We demonstrate that Picky analysis is an effective tool for comprehensive detection of SVs in cancer genomes from long-read data.

Transcriptional responses of the marine green alga Micromonas pusilla and an infecting prasinovirus under different phosphate conditions.

Prasinophytes are widespread marine algae for which responses to nutrient limitation and viral infection are not well understood. We studied the picoprasinophyte, Micromonas pusilla, grown under phosphate-replete (0.65 ± 0.07 d-1 ) and 10-fold lower (low)-phosphate (0.11 ± 0.04 d-1 ) conditions, and infected by the phycodnavirus MpV-SP1. Expression of 17% of Micromonas genes in uninfected cells differed by >1.5-fold (q < 0.01) between nutrient conditions, with genes for P-metabolism and the uniquely-enriched Sel1-like repeat (SLR) family having higher relative transcript abundances, while phospholipid-synthesis genes were lower in low-P than P-replete. Approximately 70% (P-replete) and 30% (low-P) of cells were lysed 24 h post-infection, and expression of ≤5.8% of host genes changed relative to uninfected treatments. Host genes for CAZymes and glycolysis were activated by infection, supporting importance in viral production, which was significantly lower in slower growing (low-P) hosts. All MpV-SP1 genes were expressed, and our analyses suggest responses to differing host-phosphate backgrounds involve few viral genes, while the temporal program of infection involves many more, and is largely independent of host-phosphate background. Our study (i) identifies genes previously unassociated with nutrient acclimation or viral infection, (ii) provides insights into the temporal program of prasinovirus gene expression by hosts and (iii) establishes cell biological aspects of an ecologically important host-prasinovirus system that differ from other marine algal-virus systems.

Culture-independent analysis of liver abscess using nanopore sequencing.

The identification of microbial species has depended predominantly upon culture-based techniques. However, the difficulty with which types of organisms are cultured implies that the grown species may be overrepresented by both cultivation and plate counts. In recent years, culture-independent analysis using high-throughput sequencing has been advocated for use as a point-of-care diagnostic tool. Although it offers a rapid and unbiased survey to characterize the pathogens in clinical specimens, its accuracy is reduced by the high level of contamination of human DNA. In this paper, we propose using a culture-independent analysis for a Klebsiella pneumoniae clinical strain within a liver abscess using nanopore sequencing. Owing to the highly-contaminated cell population within a liver abscess, we managed to reduce the confounding effects of human DNA through the use of DNase and differential centrifugation. Genomic DNA was sequenced through the use of Nanopore MinION sequencer and analyzed using a suite of bioinformatics approaches. K. pneumoniae was successfully identified along with antibiotic-resistant genes. Our results indicate that, by integrating real-time nanopore sequencing and bioinformatics software, real-time pathogen identification in a liver abscess can be achieved.

A statistical method for detecting differentially expressed SNVs based on next-generation RNA-seq data.

In this article, we propose a new statistical method-MutRSeq-for detecting differentially expressed single nucleotide variants (SNVs) based on RNA-seq data. Specifically, we focus on nonsynonymous mutations and employ a hierarchical likelihood approach to jointly model observed mutation events as well as read count measurements from RNA-seq experiments. We then introduce a likelihood ratio-based test statistic, which detects changes not only in overall expression levels, but also in allele-specific expression patterns. In addition, this method can jointly test multiple mutations in one gene/pathway. The simulation studies suggest that the proposed method achieves better power than a few competitors under a range of different settings. In the end, we apply this method to a breast cancer data set and identify genes with nonsynonymous mutations differentially expressed between the triple negative breast cancer tumors and other subtypes of breast cancer tumors.

Identifying Aspects of the Post-Transcriptional Program Governing the Proteome of the Green Alga Micromonas pusilla.

Micromonas is a unicellular motile alga within the Prasinophyceae, a green algal group that is related to land plants. This picoeukaryote (<2 µm diameter) is widespread in the marine environment but is not well understood at the cellular level. Here, we examine shifts in mRNA and protein expression over the course of the day-night cycle using triplicated mid-exponential, nutrient replete cultures of Micromonas pusilla CCMP1545. Samples were collected at key transition points during the diel cycle for evaluation using high-throughput LC-MS proteomics. In conjunction, matched mRNA samples from the same time points were sequenced using pair-ended directional Illumina RNA-Seq to investigate the dynamics and relationship between the mRNA and protein expression programs of M. pusilla. Similar to a prior study of the marine cyanobacterium Prochlorococcus, we found significant divergence in the mRNA and proteomics expression dynamics in response to the light:dark cycle. Additionally, expressional responses of genes and the proteins they encoded could also be variable within the same metabolic pathway, such as we observed in the oxygenic photosynthesis pathway. A regression framework was used to predict protein levels from both mRNA expression and gene-specific sequence-based features. Several features in the genome sequence were found to influence protein abundance including codon usage as well as 3' UTR length and structure. Collectively, our studies provide insights into the regulation of the proteome over a diel cycle as well as the relationships between transcriptional and translational programs in the widespread marine green alga Micromonas.

Integrative epigenomic analysis reveals unique epigenetic signatures involved in unipotency of mouse female germline stem cells.

BACKGROUND: Germline stem cells play an essential role in establishing the fertility of an organism. Although extensively characterized, the regulatory mechanisms that govern the fundamental properties of mammalian female germline stem cells remain poorly understood. RESULTS: We generate genome-wide profiles of the histone modifications H3K4me1, H3K27ac, H3K4me3, and H3K27me3, DNA methylation, and RNA polymerase II occupancy and perform transcriptome analysis in mouse female germline stem cells. Comparison of enhancer regions between embryonic stem cells and female germline stem cells identifies the lineage-specific enhancers involved in germline stem cell features. Additionally, our results indicate that DNA methylation primarily contributes to female germline stem cell unipotency by suppressing the somatic program and is potentially involved in maintenance of sexual identity when compared with male germline stem cells. Moreover, we demonstrate down-regulation of Prmt5 triggers differentiation and thus uncover a role for Prmt5 in maintaining the undifferentiated status of female germline stem cells. CONCLUSIONS: The genome-wide epigenetic signatures and the transcription regulators identified here provide an invaluable resource for understanding the fundamental features of mouse female germline stem cells.

Evidence-based green algal genomics reveals marine diversity and ancestral characteristics of land plants.

BACKGROUND: Prasinophytes are widespread marine green algae that are related to plants. Cellular abundance of the prasinophyte Micromonas has reportedly increased in the Arctic due to climate-induced changes. Thus, studies of these unicellular eukaryotes are important for marine ecology and for understanding Viridiplantae evolution and diversification. RESULTS: We generated evidence-based Micromonas gene models using proteomics and RNA-Seq to improve prasinophyte genomic resources. First, sequences of four chromosomes in the 22 Mb Micromonas pusilla (CCMP1545) genome were finished. Comparison with the finished 21 Mb genome of Micromonas commoda (RCC299; named herein) shows they share ≤8,141 of ~10,000 protein-encoding genes, depending on the analysis method. Unlike RCC299 and other sequenced eukaryotes, CCMP1545 has two abundant repetitive intron types and a high percent (26 %) GC splice donors. Micromonas has more genus-specific protein families (19 %) than other genome sequenced prasinophytes (11 %). Comparative analyses using predicted proteomes from other prasinophytes reveal proteins likely related to scale formation and ancestral photosynthesis. Our studies also indicate that peptidoglycan (PG) biosynthesis enzymes have been lost in multiple independent events in select prasinophytes and plants. However, CCMP1545, polar Micromonas CCMP2099 and prasinophytes from other classes retain the entire PG pathway, like moss and glaucophyte algae. Surprisingly, multiple vascular plants also have the PG pathway, except the Penicillin-Binding Protein, and share a unique bi-domain protein potentially associated with the pathway. Alongside Micromonas experiments using antibiotics that halt bacterial PG biosynthesis, the findings highlight unrecognized phylogenetic complexity in PG-pathway retention and implicate a role in chloroplast structure or division in several extant Viridiplantae lineages. CONCLUSIONS: Extensive differences in gene loss and architecture between related prasinophytes underscore their divergence. PG biosynthesis genes from the cyanobacterial endosymbiont that became the plastid, have been selectively retained in multiple plants and algae, implying a biological function. Our studies provide robust genomic resources for emerging model algae, advancing knowledge of marine phytoplankton and plant evolution.

Molecular portraits of epithelial, mesenchymal, and hybrid States in lung adenocarcinoma and their relevance to survival.

Epithelial-to-mesenchymal transition (EMT) is a key process associated with tumor progression and metastasis. To define molecular features associated with EMT states, we undertook an integrative approach combining mRNA, miRNA, DNA methylation, and proteomic profiles of 38 cell populations representative of the genomic heterogeneity in lung adenocarcinoma. The resulting data were integrated with functional profiles consisting of cell invasiveness, adhesion, and motility. A subset of cell lines that were readily defined as epithelial or mesenchymal based on their morphology and E-cadherin and vimentin expression elicited distinctive molecular signatures. Other cell populations displayed intermediate/hybrid states of EMT, with mixed epithelial and mesenchymal characteristics. A dominant proteomic feature of aggressive hybrid cell lines was upregulation of cytoskeletal and actin-binding proteins, a signature shared with mesenchymal cell lines. Cytoskeletal reorganization preceded loss of E-cadherin in epithelial cells in which EMT was induced by TGFß. A set of transcripts corresponding to the mesenchymal protein signature enriched in cytoskeletal proteins was found to be predictive of survival in independent datasets of lung adenocarcinomas. Our findings point to an association between cytoskeletal and actin-binding proteins, a mesenchymal or hybrid EMT phenotype and invasive properties of lung adenocarcinomas.

Lineage-specific chromatin signatures reveal a regulator of lipid metabolism in microalgae.

Alga-derived lipids represent an attractive potential source of biofuels. However, lipid accumulation in algae is a stress response tightly coupled to growth arrest, thereby imposing a major limitation on productivity. To identify transcriptional regulators of lipid accumulation, we performed an integrative chromatin signature and transcriptomic analysis to decipher the regulation of lipid biosynthesis in the alga Chlamydomonas reinhardtii. Genome-wide histone modification profiling revealed remarkable differences in functional chromatin states between the algae and higher eukaryotes and uncovered regulatory components at the core of lipid accumulation pathways. We identified the transcription factor, PSR1, as a pivotal switch that triggers cytosolic lipid accumulation. Dissection of the PSR1-induced lipid profiles corroborates its role in coordinating multiple lipid-inducing stress responses. The comprehensive maps of functional chromatin signatures in a major clade of eukaryotic life and the discovery of a transcriptional regulator of algal lipid metabolism will facilitate targeted engineering strategies to mediate high lipid production in microalgae.

Marine algae and land plants share conserved phytochrome signaling systems.

Phytochrome photosensors control a vast gene network in streptophyte plants, acting as master regulators of diverse growth and developmental processes throughout the life cycle. In contrast with their absence in known chlorophyte algal genomes and most sequenced prasinophyte algal genomes, a phytochrome is found in Micromonas pusilla, a widely distributed marine picoprasinophyte (<2 µm cell diameter). Together with phytochromes identified from other prasinophyte lineages, we establish that prasinophyte and streptophyte phytochromes share core light-input and signaling-output domain architectures except for the loss of C-terminal response regulator receiver domains in the streptophyte phytochrome lineage. Phylogenetic reconstructions robustly support the presence of phytochrome in the common progenitor of green algae and land plants. These analyses reveal a monophyletic clade containing streptophyte, prasinophyte, cryptophyte, and glaucophyte phytochromes implying an origin in the eukaryotic ancestor of the Archaeplastida. Transcriptomic measurements reveal diurnal regulation of phytochrome and bilin chromophore biosynthetic genes in Micromonas. Expression of these genes precedes both light-mediated phytochrome redistribution from the cytoplasm to the nucleus and increased expression of photosynthesis-associated genes. Prasinophyte phytochromes perceive wavelengths of light transmitted farther through seawater than the red/far-red light sensed by land plant phytochromes. Prasinophyte phytochromes also retain light-regulated histidine kinase activity lost in the streptophyte phytochrome lineage. Our studies demonstrate that light-mediated nuclear translocation of phytochrome predates the emergence of land plants and likely represents a widespread signaling mechanism in unicellular algae.

A search for novel cancer/testis antigens in lung cancer identifies VCX/Y genes, expanding the repertoire of potential immunotherapeutic targets.

Cancer/testis (CT) antigens are potential immunotherapeutic targets in cancer. However, the expression of particular antigens is limited to a subset of tumors of a given type. Thus, there is a need to identify antigens with complementary expression patterns for effective therapeutic intervention. In this study, we searched for genes that were distinctly expressed at a higher level in lung tumor tissue and the testes compared with other nontumor tissues and identified members of the VCX/Y gene family as novel CT antigens. VCX3A, a member of the VCX/Y gene family, was expressed at the protein level in approximately 20% of lung adenocarcinomas and 35% of squamous cell carcinomas, but not expressed in normal lung tissues. Among CT antigens with concordant mRNA and protein expression levels, four CT antigens, XAGE1, VCX, IL13RA2, and SYCE1, were expressed, alone or in combination, in about 80% of lung adenocarcinoma tumors. The CT antigen VCX/Y gene family broadens the spectrum of CT antigens expressed in lung adenocarcinomas for clinical applications.